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The stability of a vortex glass phase with quasi-long-range positional order is examined for a 
disordered layered superconductor. The role of topological defects is investigated using a detailed 
scaling argument, supplemented by a variational calculation. The results indicate that topological 
order is preserved for a wide range of parameters in the vortex glass phase. The extent of the 
stability regime is given in terms of a simple Lindemann-like criterion. 
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It is well-known that the Abrikosov flux lattice in a 
type-II superconductor is unstable to point disorder be- 
yond the Larkin length The nature of the flux ar- 
ray at larger scales has been a subject of intense stud- 
ies It has been conjectured that the flux array is 
collectively pinned, forming a vortex glass (VG) phase 
1^-^ with zero linear resistivity at low temperatures. 
This conjecture is supported by a number of experiments 
on disordered samples of high-Tc superconductors [[/[-p^ , 
where a continuous transition to a phase with zero lin- 
ear resistivity was found upon cooling. On the other 
hand. Bitter-decoration jl^, neutron scattering [Q, and 
/iSR experiments on weakly-disordered samples have 
all indicated some long-range order of the flux array, a 
characteristic usually incompatible with a glass. A com- 
mon interpretation for the observation of a flux lattice 
is a crossover effect due to the large Larkin lengths in 
weakly disordered samples. In this article, we investigate 
a different possibility, that the flux array may maintain 
its positional long-range order much beyond the Larkin 
length, in spite of pinning by point disorders. 

Such a possibility is indeed realized in a model of 
dislocation-free flux line array in random media |p|,p^. 
This model is very similar to the randomly-pinned 
charge-density waves and the random-field XY model 
which have been studied extensively in the past decades 
iP^-pTf. A variety of approximate methods have been 
used to obtain the conclusion that point disorders lead 
to a glass phase with only logarithmic fluctuations in the 
transverse displacement of the flux array. This implies 
the existence of quasi-long-range positional order in the 
glass phase of the dislocation- free flux array Q . 

Recently, Giamarchi and LeDoussal suggested QjlSj 
that such a topologically-ordered glass may actually ex- 
ist as a stable thermodynamic phase for some range of 
parameters in the cuprate superconductors. A related 
numerical study of the random field-XY model by Gin- 



ever, the issue of spontaneous formation of topological 
defects (i.e., dislocation loops) involves complicated in- 
terplay between elasticity and disorders, and has so far 
not been addressed quantitatively. In this article, we in- 
vestigate this issue by studying a model of flux lines con- 
flned in the planes of a layered superconductor [Fig. 1]. 
Our model allows for the formation of dislocation loops 
and is amenable to analytic studies. We first present 
a detailed scaling argument, which yields suppression of 
large dislocation loops at finite fugacities. This result 
is then supplemented by a variational calculation, from 
which we obtain a Lindemann-like criterion giving the 
size of the stability regime for the topologically-ordered 
VG. Similar behaviors are expected in the usual experi- 
mental situation of fiux lines perpendicular to the layers. 




gras and Huse |19| further supported this scenario. How- 



FIG 1: Flux line array in a layered superconductor in a 
parallel magnetic field. 

We consider a strongly layered impure superconductor 
in a parallel magnetic field [see Fig. 1]. The supercon- 
ducting layers provide a sufficiently strong confining po- 
tential for the (Josephson-like) vortex lines which exist in 
the interlayer spacing. We shall exclude the possibility of 
the lines crossing the superconducting layers. [For fields 
parallel to the ab-planes of the Bi-compound, typical vor- 
tex kink energies are of the order 10^(1 — T/T^) K] This 
amounts to limiting the vortex displacement field from 
two components in an isotropic sample to one component 
(i.e., parallel to the layers). For simplicity we focus our 
study mostly on the dilute limit where the inter-vortex 



1 



spacings l± (inter-layer), Z|| (intra-layer) exceed the mag- 
netic penetration depths Xab, Ac respectively. The impli- 
cations of our results on the dense limit are straightfor- 
ward and will be discussed below. 

A well-established analytic description for a single 
layer of vortex lines (for length scales exceeding is 
given by the Hamiltonian |^ 

f3n2D[cPj,W,]= 1^ |^(V||(/.,f-M/,[</.,(r||),r||]|, 

(1) 

where ^j(r||) describes the in-plane displacement of the 
vortex lines in the j^^ layer and K is an (isotropized) in- 
plane elastic constant. Pinning effects due to point disor- 
der are described by the random potential [(/)j(r||), ry], 
with the second moment 

W^,[0„r||]T4^,[0;.,r||] =.go'cos[(/)-0']<52(r|| -r||), (2) 

where the overbar denotes disorder average, go character- 
izes the (bare) strength of the random potentials, and the 
cosine captures the discrete nature of the vortex lines . 
With many layers stacked next to each other, the Hamil- 
tonian for the whole system is 

= E I /3^2D + / V,[cl,,+,m)- </.,(r||)] I , (3) 



where WjWj' = Sjj'WjWj since the bare random poten- 
tials in different layers are uncorrelated, and the interac- 
tion has the form 

1^,M = -Mcos[0]. (4) 

The expression (^) can be regarded as the repulsive mag- 
netic interaction energy between the lowest harmonics of 
density fluctuations between vortex lines in "adjacent" 
layers, a valid approximation in the dilute limit ]20| ]. The 
coupling constant is related to the shear modulus of the 
flux line lattice. The main feature of (^) is that it goes 
beyond the elastic approximation, as it allows for dislo- 
cation loops between adjacent layers [see Fig. 2]. 




FIG 2: Elastic rearrangement of each planar flux array 
can be represented by a number of vortex loops. An 
example is shown on the left (shaded region). Aligned 
vortex loops in successive layers form a vortex sheet as 
shown on the right. The boundary of the sheet (dark 
loop on top) is a dislocation loop. 

In what follows, we first study the phase diagram of 
the system in terms of the parameters ^ and K. If the 
vortex layers are uncoupled, i.e., /it = 0, then each layer 
undergoes separately a glass transition at a critical value 
Kc — l/47r Pjl^, with a vanishing linear resistivity in 
the low temperature phase {K > Kc) and ohmic behav- 
ior in the high temperature phase {K < Kc) j2|j2^. The 
physical range of interest corresponds to K ^ Kc. When 
coupling the layers together via Eq. (^), a competition 
arises between the tendency of vortex lines in each layer 
to lower their free energy by adapting to the disorder in 
their own layer, against the tendency to minimize the in- 
terlayer coupling by bringing neighboring vortex layers 
into registry. The limit of very weak coupling (fi ^ 1) 
may be studied using perturbation theory; it is straight- 
forward to find that weak coupling is irrelevant at large 
scales. 

In the limit of very large coupling /i — > cx), the inter- 
action potential V^[0] may be replaced by the quadratic 
form, §0^, which describes an elastic (i.e., dislocation- 
free) coupling in the direction perpendicular to the lay- 
ers. This is just the anisotropic (and one-component) ver- 
sion of the VG considered previously in Refs. P, p^ , p^ . 
After introducing a continuum description in terms of 
r = (r||,r_L) and rescahng r± = {j ■ 1\^^/JjI\JK , we get 
an isotropic 3D elastic Hamiltonian 

/?H3D =yd'r{|(V0)2-W^[0(r),r]} (5) 

with an effective elastic constant 7 = ^/JlK and a random 
potential iy[(/)(r),r] with iy[0, r]M^[0, 0] = g"^ cos[cj)]5^ {v) , 
where g^ — g^^Jfl/K. From various methods including 
position-space RG |jl^ , functional RG , and the vari- 
ational Ansatz with replica-symmetry breaking |]l^ , |l7| , 
one finds that the system (||) forms a glass phase with 

mv)-<j>{v')Y)^^ = 2A\og{\v-v'\/t) (6) 

beyond the positional correlation length £ ~ 7^/5^, with 
A being a universal number of 0(1). [In this model, £ is 
equivalent to the Larkin length, although in more realis- 
tic models including higher harmonics of density fluctu- 
ations, £ exceeds the Larkin length [^.]The logarithmic 
fluctuation in displacement leads to quasi-long-range or- 
der and an (algebraic) Bragg peak at reciprocal lattice 
vector 27r//|| . This phase has since been referred to as the 
"Bragg glass" [Q. In 3D, there is a large elastic energy 
cost, of the order (7/2) J d^r{{'V4>)'^)^Y) ~ 7-^' ^'^^ loga- 
rithmic fluctuations in a volume of the order L^. This en- 
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ergy is compensated by the disorder energy gained from 
the anomalous displacement of the flux array. Thus 

AS - 7L (7) 

gives the sample-to-sample free energy fluctuation of the 
Bragg glass. Also, due to the existence of a statistical 
symmetry |^,^ , it is known that the disorder-averaged 
responses of the system to various elastic deformations 
are identical to those of the pure elastic system, i.e., 
Eq. (||) with W = 0. In particular, there is a non-zero 
response to shear. 

Given the above properties of the Bragg glass, which 
exists so far only in the unphysical limit /i — > cxd, our first 
task is to determine whether it can persist at a finite /i, 
i.e., whether the system is stable to the spontaneous for- 
mation of dislocation loops on length scales much larger 
than the correlation length. To investigate this possibil- 
ity, we divide the system into two halves (within which 
the layers are elastically coupled), and allow dislocation 
loops to form in the contact plane, say between the jj'' 
and (jo + 1)*'' layers. Analytically, this is implemented 
by using the following interaction energy in Eq. (^: 

V^ict:] ^^Cj,' if J^JO 

= -fi' cos[(j)] if j^jo, (8) 

where fj.' fa fi approximates (^). It is instructive to con- 
sider arbitrary values of ^' . Let us compare the free en- 
ergy difference between the completely coupled and com- 
pletely decoupled limit. If the two halves of the system 
are decoupled, then each forms a Bragg glass, and the 
configuration of the flux array in each half is individually 
optimized. But if the two halves are tied together, then 
the constraint across the contact plane forces a global re- 
optimization of the flux array, resulting in a higher free 
energy for each half. The typical free energy increase 
due to the constraint is given by the sample-to-sample 
free energy fluctuation of the Bragg glass, AE ^ jL. 
The expression (0) is used here because each half of the 
system must be completely re-optimized given the con- 
straint, as if the conflguration of the random potential 
has been completely changed. The difference in the opti- 
mal configuration of the fiux array in each half resulting 
from the constraint at jo can be described by a collec- 
tion of "vortex sheets" such as the one depicted in Fig. 2. 
The boundary of a vortex sheet is a dislocation loop; it 
describes phase mismatches across the contact plane. 

To find whether or not the two half systems actually 
couple for a given /i', it is necessary to balance the cost in 
disorder energy due to coupling with the reduction in in- 
teraction energy due to phase matching. The latter can 
be readily computed for small /i' using (^). One finds 
/i'/i2 d^iX(cos[0jo+i - 0jo])3j5 fi'i^L^-^ < L since 
A > 1 Ej]. As the disorder energy cost to coupling, 
AE ~ jL, exceeds the interaction energy to be gained 
at small /j,' and large L, the two half systems will remain 



decoupled. In the large (but finite) n' limit of interest, 
the perturbative result is no longer valid. Since the en- 
ergy cost of phase mismatch is large there, we consider 
the stability of a single optimally-configured dislocation 
loop of extent L ^ £ at the contact plane of the two half 
systems that arc otherwise elastically coupled [Fig. 2]. 
The energy cost of the core of the dislocation loop due to 
the inter-layer interaction is extensive. For a stretched 
circular loop of linear size L, we expect i?coro ~ ^j.'£L. 
Here, £ appears as the "thickness" of the loop because 
the flux array is elastically coupled at smaller scales at 
low temperatures |26|. More generally, if we allow the 
dislocation loop to take on fractal shapes, say with the 
total loop length scaling as for L :p £ {D > 1 be- 
ing the fractal dimension) , then the core energy becomes 
-E-corc ~ fj! £^~^ . On the other hand, the disorder 
energy gained from the formation of a dislocation loop 
of size L is just the energy gained from the formation 
of a vortex sheet Esheet{L), resulting from the elastic 
deformation of the half-system [see Fig. 2]. Esheet{L) 
clearly cannot exceed AE{L) which is the disorder en- 
ergy gained from complete optimal elastic rearrangement 
of the half-system at scale L. Assuming scaling of vortex 
sheet energies, Esheet{L) ^ it follows that lo <\. 

The actual value of the exponent to depends on the 
structure of the dislocation loop we allow, i.e., on the 
fractal dimension D. We expect that the upper bound 
(w = 1) may only be reached if the structure of the as- 
sociated vortex sheet is similar to those occurring natu- 
rally in complete elastic rearrangement of the half system. 
The latter can be deduced as follows: Denote the dif- 
ference in the configuration before/after rearrangement 
by ((9(r). The vortex sheets are then the equal-(^ con- 
tours of (/'(r), and the associated dislocation loops are 
the contours of tp(Y\\,r^ = jol±)- The relationship be- 
tween a rough "landscape" and the fractal geometry of 
its contours have recently been examined [ p7[ . For a 
logarithmically-rough landscape (p [see Eq. (o)], an ex- 
act calculation yields D = 3/2 [pTj-pOf. Thus, we expect 
Lu <1 for D = 3/2, and w < 1 for D < 3/2. The total 
energy of the dislocation loop 

Sloop = Score (i) - S,heet(i) ^ /i'^'-^L^ - t'^"^""^ 

(9) 

does not admit a stable solution with L 3> ^ for large 
/i' « fi. Hence the Bragg glass is stable to the spon- 
taneous formation and proliferation of large dislocation 
loops. This conclusion is further supported by a sys- 
tematic renormalization-group analysis, the details of 
which will be given elsewhere The possibility of a 

marginally stable Bragg glass for weakly-disorder sample 
was first suggested in Ref. [Q, based on the assertions 
that AE ^ gL and Score cL, where g is the bare disor- 
der strength which can be made arbitrary small and c is a 
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given number . The above analysis indicates that the 
dislocation loops are much more strongly suppressed at 
low temperatures by the anomalously large core energy. 

Next, we investigate the extent of the stability regime 
for the Bragg glass. As already shown, the Bragg glass 
cannot exist in the limit of very weak inter-layer cou- 
pling. We expect the maximum possible extent of the 
stability regime to be the point where the disorder is so 
strong that the correlation length £ becomes of the or- 
der of the average vortex spacing a. Beyond this point, 
a topologically-ordered flux array cannot even form at 
the smallest scale, and the entire collective pinning pic- 
ture breaks down. To find the actual size of the stabil- 
ity regime for the layered system (^), we consider quasi 
two-dimensional in-plane fluctuations i.e., on the shortest 
scale in the direction perpendicular to the layers. Analyt- 
ically, we apply a variational treatment, with the varia- 
tional Hamiltonian obtained by replacing the interaction 

potential in (^) by the quadratic form — ^(jj"^. 
V describes an elastic (i.e., dislocation- free) coupling in 
the direction perpendicular to the layers. The param- 
eter Jl has the meaning of an effective shear modulus 
and may be determined self-consistently within the vari- 
ational treatment. The minimization of the variational 
free energy with respect to Jl yields the self-consistency 
equation 

Jl = /i(cos[(?!)j+i(r||) - (t>j{^)])~. (10) 

This is evaluated using a Gaussian approximation which 
can be justified |3^] in a controlled fashion and should 
be reasonably accurate. {<p(f))~ contains contributions 

from (i) the quasi 2D VG regime {(j)4>)2D which domi- 
nates for Jl fa 0; (ii) the 3D VG regime {4'(j})^ij; and (iii) 
thermal fluctuations on scales smaller than the correla- 
tion length £ for large Jl. Using Eq. (^ for and 
using (0^)^= 2(1 -I- Q!)log(L) (where a fa [l\\go/KY 
from Refs. p^ , p7t ), the following results are obtained 
The self-consistency equation has a stable solution with 
nonzero shear modulus only for |^ /i > /^c = c^K/£'^. 
For n < lie, the system "melts" into a stack of decoupled 
2D VG's, distinguished from the Bragg glass by a vanish- 
ing shear modulus Jl — fs^ . The transition at fi = fic 
is first order, with a jump in Jl |3C| ]. Our variational 
calculation yields a prefactor c ~ 5 which depends very 
weakly on temperatures, as long as we are away from the 
melting temperature of the pure system [^o| . 

It is illustrative to express fi and K in terms of the 
correlation length [of the anisotropic system (||)] in the 
_L-direction, £± = ijlI\/K£. The above stability condi- 
tion becomes 

£a^>c-Ia_. (11) 

Thus, we find the Bragg glass to be very stable, reaching 
within a factor of 10 of the maximum extent of stability 



(i.e., £ K. I) beyond which the very concept of collective 
pinning breaks down! As variational calculations tend to 
under-estimate fluctuations, we expect the actual value 
of c to be somewhat larger. Nevertheless, the condition 
( prj ) appears to be quite general, and may be viewed as 
the disordered-analog of the Lindemann criterion. 

Applying this Lindcmann-like criterion, we can expect 
the following low-temperature behaviors for the usual ex- 
perimental situation of flux lines perpendicular to the 
CuO planes; a detailed account including the relevance 
to recent experiments will be given elsewhere (a) 
In the dilute limit, the Bragg glass phase is stable if the 
transverse correlation length £{B) exceeds (a few times) 
the inter- vortex spacing 1{B). (b) The same conclusion 
is obtained in the dense limit [A if £ > A. (c) In the 
opposite limit of ^ < £ ^ A ^ oo, the long-ranged mag- 
netic interaction leads to a much stronger glass phase, 
with Ai? ~ [0, if the dislocation loops are forbid- 
den. This however provides a much larger disorder en- 
ergy to be gained from the formation of dislocation loops 
[i.e.,Eq. (^ with w = 2], leading to the proliferation of 
large loops. This regime may belong to the universal- 
ity class of the (long-ranged) Gauge glass |^^, or it may 
simply be a viscous line liquid, (d) For the physical case 
of finite A's, dislocation loops can form at scales A and 
below, and we expect the system to remain topologically- 
ordered at large scales. 

In summary, we find strong evidence supporting the 
existence of a topologically-ordered vortex glass for the 
layered superconductors. The stability regime is surpris- 
ingly large and is given by a Lindemann-like criterion. 
Beyond this regime, the system is dominated by strong 
disorders^ and its properties are not known. It may sim- 
ply melt into a viscous line liquid, or it may form another 
type of glass phase | |T^ . 
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